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Abstract 

A great deal of progress has been made in the 
modeling of aerodynamically generated sound for ro- 
tors over the past decade. The Ffowcs Williams- 
Hawkings (FW-H ) equation has been the founda- 
tion for much of the development. Both subsonic 
and supersonic quadrupole noise formulations have 
been developed for the prediction of high-speed im- 
pulsive noise. In an effort to eliminate the need to 
compute the quadrupole contribution, the FW-H 
has also been utilized on permeable surfaces sur- 
rounding all physical noise sources. Comparison of 
the Kirchhoff formulation for moving surfaces with 
the FW-H equation have shown that the Kirchhoff 
formulation for moving surfaces can give erroneous 
results for aeroacoustic problems. 

1. Introduction 

Noise has been an undesirable byproduct of 
aerospace vehicles from the time of early aircraft un- 
til now. Originally, aircraft noise was not much of 
an issue because of the overarching requirement of 
improving vehicle performance. As aerospace tech- 
nology has matured, more resources have been de- 
voted to the reduction of aerodynamically generated 
sound. With the maturation of the aerospace tech- 
nology already achieved, both the public and regu- 
latory bodies have focused their concern on safety, 
emissions, and noise rather than either performance 
or efficiency. This situation is really a credit to 
the success of past generations of aerospace design- 
ers and engineers. The challenge now facing us is 
that anticipated increases in the utilization of air 
travel will also bring unacceptable increases in air- 
craft noise — if nothing is done. Fortunately, a great 
deal of progress has been made both in understand- 
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ing the noise-generation mechanisms and in devel- 
oping first-principles models for prediction of the 
sound. 

Rotorcraft are inherently complex aeromechanical 
vehicles, and hence have lagged fixed-wing aircraft 
in the understanding of the mechanisms responsible 
for rotor noise and accurate and efficient prediction 
methods based upon the fluid physics (as opposed 
to empirical methods). A number of high qual- 
ity experimental tests, concentrated national pro- 
grams and rigorous theoretical developments have 
greatly expanded both our understanding of rotor 
noise sources and our ability to predict rotor noise. 
The maturation of computational fluid dynamics 
(CFD) and its application to rotor aerodynamics ha s 
been an enabling step necessary for accurate, first- 
principles noise prediction. Nevertheless, rotor aero- 
dynamics is often one of the weak links leading to 
unsatisfactory acoustic computations. 

In this paper, the focus will be on the aerody- 
namically generated sound of rotors and the recent 
advances in rotor noise prediction. In 1994, Brent- 
ner and Farassat published a status report on heli- 
copter noise prediction. Their report provides a his- 
torical perspective and method assessment that is a 
useful backdrop for this paper. A significant body 
of work has been completed since 1994, therefore, 
it is the intent of this paper to put both theoreti- 
cal and computational developments in perspective. 
This paper places some emphasis on theoretical de- 
velopment and interpretation because these always 
help explain the reason various computational meth- 
ods were chosen. Advances in source noise prediction 
are the primary focus of this paper even though sys- 
tem noise prediction is recognized as an important 
area that still lacks sufficient capabilities. 

2. Theoretical Background: FW-H equation 

The problem of aerodynamically generated sound 
is governed by the conservation laws of mass, mo- 
mentum and energy. The FW-H equation 1 is an 
exact rearrangement of the continuity equation and 
the Navier-Stokes equations into the form of an in- 
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homogeneous wave equation with two surface source 
terms and a volume source term. The FW-H equa- 
tion is the most general form of the Lighthill acous- 
tic analogy 2 * and is appropriate for predicting the 
noise generated by the complex motion of helicopter 
rotors. Today almost all deterministic rotor noise 
predictions are based on time-domain integral for- 
mulations of the FW-H equation. In this section, 
the FW-H equation will be examined to provide the 
background necessary to understand the recent ad- 
vances in rotor noise prediction. The close relation- 
ship between the FW-H equation and the Kirchhoff 
formulation for moving surfaces will be discussed 
later. 


equation is valid for the entire space — both inside 
and outside of the body. The generalized momentum 
equation can be written 

dpui BpuiUj dPij _ 

dt dxj dxj 

{PUi-r £ + {PUiUj + p ij)^) 5 U) ■ ( 5 ) 

Now by taking the time derivative of equation (4) 
and subtracting the divergence of equation (5), fol- 
lowed with some rearranging, the FW-H equation 
may be written as the following inhomogeneous wave 
equation: 


2.1 Derivation 

The FW-H equation may be derived by em- 
bedding the exterior flow problem in unbounded 
space by using generalized functions to describe the 
flow field. To do this, consider a moving surface 
/(x,£) = 0 with a stationary fluid outside. The sur- 
face / = 0 is defined such that V/ = n, where n 
is a unit normal vector pointing toward the exte- 
rior of the surface. Inside / = 0 the generalized flow 
variables are defined to have their freestream values, 
i.e., 


p = p'H (/) + p 0 

( 1 ) 

pvii = puiH(f) 

( 2 ) 

= PljHW+pJij 

( 3 ) 


where the tilde indicates that the variable is a gen- 
eralized function defined throughout all space. On 
the right hand side p, piq, and Pij are the density, 
momentum, and compressive stress tensor, respec- 
tively. Freestream quantities are indicated by the 
subscript o, primed quantities represent the differ- 
ence from the freestream value (e.g., p' = p - p 0 ), 
5{j is the Kronecker delta and H(f) is the Heaviside 
function. (Note: for an inviscid fluid, P \j = pSij.) 

Substituting the generalized density and momen- 
tum definitions (l)-(2) into the continuity equation 
and collecting the derivatives of the Heaviside func- 
tion on the right-hand side yields a generalized con- 
tinuity equation written as 


dp dpui 
dt dxi 




( 4 ) 


Here the bar over the derivative operators indicate 
that generalized differentiation (i.e., differentiation 
of generalized functions) is implied. Also note that 
df/dt = -u n , df/dxi = hi and dH(f)/df = 5(f) is 
the Dirac delta function. This generalized continuity 


- - ««))<*(/)] 

+ | ] - t \{PoV n + p(u n - V n ))6(f)] ( 6 ) 

where T{j is the Lighthill stess tensor, u n is the fluid 
velocity in the direction normal to the surface / = 0 
and v n is the surface velocity in the direction normal 
to the surface. On the left hand side we use the 
notation D 2 = [(l/c 2 )(d 2 /dt 2 )\ — V 2 . 

If the fictitious surface / = 0 coincides with a solid 
surface, then the normal velocity of the fluid is the 
same as the normal velocity of the surface ( u n = v n ). 
Then equation (6) can be written 

n v<M) = 

- [{P'anAf)] + J t [(PoV n )6(f)] 

( 7 ) 

where on the left hand side we use the customary 
notation p r = c 2 p f because the observer location is 
outside of the source region. 

2.2 Interpretation 

The three source terms on the right-hand side of 
equation (7) are known as the thickness, loading, 
and quadrupole source terms, respectively. A for- 
mal solution may be obtained by using the free- 
space Green’s function and utilizing the fact that 
the FW-H equation is valid in the entire unbounded 
space. 

Note in equation (7) that the thickness and load- 
ing source terms are surface distributions of sources 
(indicated by the presence of the Dirac delta func- 
tion 5(f)). The thickness source accounts for the 
noise generated by the displacement of fluid as the 
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body passes while the loading source term accounts 
for noise that results from the unsteady motion of 
the force distribution on the body surface. The 
thickness and loading source terms have been used 
for several years in rotor noise prediction because 
they account for most of the acoustic signal when 
the flow field is not transonic. Furthermore, they do 
not require knowledge of the flow field off the blade 
(although the accurate determination of the blade- 
surface pressure is still challenging). 

The quadrupole source, on the other hand, is 
a volume distribution of sources (indicated by the 
Heaviside function H(f )). The quadrupole source 
accounts for the nonlinearities due to both the lo- 
cal sound speed variation and the finite particle 
velocity near the blade. The importance of the 
quadrupole term has long been recognized, 3,4 how- 
ever, the quadrupole source has often been neglected 
in rotor noise prediction because of the computa- 
tional demands of computing the flow field with suf- 
ficient accuracy and integrating over a volume in the 
acoustic prediction. 

It is interesting to note that the form of the various 
source terms is not unique. 5 One example of this, 
known as Isom thickness noise, is that the thick- 
ness source term in equation (7) is equivalent to a 
constant pressure c?p 0 over the body in the loading 
source term. 6 Although mathematically equivalent, 
the two formulations for thickness noise have quite 
different characteristics and robustness when inte- 
grated numerically. Much of the recent work on the 
FW-H equation has been connected with efficient 
evaluation of the quadrupole source term or develop- 
ment of a more numerically suitable, but equivalent 
source description. 


solution may be written as 

4 = j J sk^mnm dydT . m 

— oo — oo 

The next stage in developing the acoustic formula- 
tion is to integrate the Dirac delta functions 5(f) 
and 5(g ), a process that requires a change of vari- 
ables. This change of variables determines the type 
of formulation. Equation (9) can be expressed as 

4 ^ < 10 > 

/= 0 

= [ [ cdTdr ( 11 ) 

J J rsm6 
"°° /= o 

= (i2) 

F —0 

with the variable transformations (r, y$) — > (<?,/), 
(y 2 , 2/s) -> (/,$), and (r,y 3 ) -> (g,F), respec- 
tively. (See reference 7 for definitions.) The 
three formulations expressed in equations (10)-(12) 
will be termed retarded-time, collapsing-sphere, and 
emission-surface formulations, respectively. Each 
type of formulation has its own physical and geo- 
metrical interpretation. For volume source terms, 
such as the quadrupole, only retarded-time and 
collapsing-sphere formulations are found. 


2.3 Integral Formulations 

A key feature of the FW-H equation is that it 
is an inhomogeneous wave equation for the exter- 
nal flow problem that has been embedded in un- 
bounded space. Hence, an integral representation of 
the solution can be readily found using the free-space 
Green’s function. The particular formulation that is 
developed results primarily from the choice of the 
change of variables needed to analytically integrate 
the Dirac delta functions. 

To illustrate, we consider the following example. 
An inhomogeneous wave equation can be written as 


n 2 <P(x,t) = Q(x,t)8(f) 


where Q(x, t) is the source strength. Equation (8) 
is typical of the various surface-source terms in the 
FW-H equation. By using the free-space Green’s 
function 5(g) /Anr , an integral representation of the 


3. Recent Advances 

Now that the ground work has been laid, it is 
time to consider some of the recent advances in rotor 
noise prediction. This presentation focuses primarily 
on formulation development and implementation at 
NASA Langley during the 1990’s. 


3.1 High-Speed Impulsive Noise 


High-speed impulsive (HSI) noise is a particu- 
larly intense and annoying noise generated by he- 
licopter rotors in high-speed forward flight. HSI 
noise is closely associated with the appearance of 
shocks and transonic flow around the advancing 
rotor blades and is accounted for by the FW-H 
quadrupole source. Farassat and Brentner 8 have 
shown that after some manipulation of the formal 
solution the noise contribution from the quadrupole 
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may be expressed as 


4^(x,() = / / 

-oo />0 

♦ *// 


-dfldr 


3T rr Ta 


dtldr (13) 


+ c 


-oo /> 0 

3T rr - Tu 


ii 

-oo />0 




where pg(x,£) is the acoustic pressure due to the 
quadrupole source. The quantity T rr is the double 
contraction of the Lighthill stress tensor with fi 
and fj ; where r* are the components of the unit vec- 
tor in the radiation direction. In addition, dfl is an 
element of the surface g = 0, which is known as the 
collapsing sphere; hence, equation (13) is known as 
the collapsing sphere formulation. Equation (13) is 
interpreted as an inner integration over the sphere 
of radius r = c(t - r) centered at the observer. The 
outer integration of -oo < r < t sums the contribu- 
tion from all spheres. 

An integration over the entire collapsing sphere 
surface is not actually necessary in the inner inte- 
grals of equation (13) because the Lighthill stress 
tensor Tij vanishes away from the source region. For 
an in-plane observer in the far field, the collapsing 
sphere can be locally approximated by a right circu- 
lar cylinder, as shown in figure 1. Because the ob- 
server is assumed to be in the rotor plane (precisely 
where HSI noise has maximum directivity) , integra- 
tion in the normal direction to the rotor plane can 
be done independent of the observer position. Yu et 
al. 3 were the first to use this far-field approximation 
for the evaluation of quadrupole noise; however, sev- 
eral additional approximations were made to both 
the quadrupole source strength and the acoustic in- 
tegrals. These additional approximations are no 
longer necessary. 

In the work of Brentner and Holland, 9 the inte- 
gration over the approximate collapsing sphere sur- 
face is carried out in two parts. First, integration in 
the direction normal to the rotor disk is performed. 
They define the quadrupole source strength on the 
rotor plane as 


Qij= J Tijdz (14) 

/> o 

where z is understood to be in the direction normal 
to the rotor disk and the 2 integration is only done 


Contours of 
source strength 



Approximate 
integration surface 




Spherical 

integration 

surface 


Figure 1. Far-field approximation of the collapsing 
sphere. Although not shown, the observer is three 
rotor radii to the right of the rotor blade. 


outside of the rotor blade; Qij is nonzero only over 
a region near the rotor blade planform, and extends 
ahead of the leading edge, behind the trailing edge, 
and off the blade tip. By using relation (14), equa- 
tion (13) may now be written as 

4^(x,f) = ~ I I ~r dTdr 

/+= 0 
9=0 

+ !/ / 

-°° /+= 0 
9—0 

+ cf J 

/+= 0 
9=0 


where / + = 0 represents the rotor-disk plane. The 
intersection of the collapsing sphere with the rotor 
plane results in a curve for which we use the notation 

r. 

3.1.1 Subsonic Formulation 

In the development of a subsonic quadrupole for- 
mulation, equation (15) is transformed from a col- 
lapsing sphere formulation to a retarded-time formu- 
lation (as was discussed previously) . When the time 
derivatives are taken inside the retarded-time inte- 
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grals, 10 the quadrupole formulation may be written 

4^(x,o= j [§^ + f? + 7rL<w (»•) 


where 

K r 1 = 


and 

K r 3 


/+=0 


Qrr M r Q rr -f 3M r Qr 


(1 - M r )3 
3M r 2 Q rr 


(1 - M r )4 


A r2 = 


(1 - M r )5 

Qit Mr M rQii 


(1 - M r ) 2 (1 - 

Mr) 2 

3[(1 - 

M 2 )Qrr - 2M r Q Mr - MiMiQrr ] 


(1 - M r Y 


6M r (l 

- M 2 )Q rr 


(1- 

-Mrf 


2Q MM 

- (1 - M 2 )Qu 

6(1 -M 2 )Q Mr 


(1 - M r Y 
3(1 - M 2 ) 2 Q rr 


(1 - M r Y 


(1 - M r f 


Equation (16), together with the definitions of K r 1} 
K r 2, and K r 3, are referred to as formulation Q1A. 
Formulation Q1A does not require numerical time 
differentiation of the integrals and, as a retarded- 
time formulation, is well suited for subsonic source 
motion. Aside from the problem geometry, only 
the time-dependent value of Qij is required as in- 
put. Brentner 10 has implemented formulation Q1A 
in a new version of the WOP WOP noise prediction 
code n now called WOPWOP+. (Full details of the 
derivation are given in reference 10.) 

To demonstrate the forward-flight capability of 
the WOPWOP+ code, Brentner 10 made a com- 
parison of predicted and measured results for a 
four-blade swept-tip rotor tested in the Duits- 
Nederslandse Windtunnel (DNW). For this compar- 
ison, a microphone located in the rotor plane at a 
rotor azimuth of ifr = 150° was utilized. The ex- 
periment is described in the report by Visintainer et 
al. 12 

The full potential solver FPRBVI 13 was used to 
compute the unsteady flow field around the rotor. 
The CFD solution was stored at every degree of rotor 
azimuth for the quadrupole source strength compu- 
tation. The results of the forward-flight noise predic- 
tion are shown in figure 2. The experimental data is 
compared with the predicted acoustic pressure; the 



Figure 2. Contemporary design, four-blade model 
rotor operating in forward flight; /x = 0.32 and 
Mat = 0.933. ■ experimental data; pre- 
dicted acoustic pressure; quadrupole compo- 

nent of predicted acoustic pressure. 


quadrupole contribution is also shown to indicate its 
relative magnitude. Although the CFD calculation 
used a rather coarse grid, the agreement is good. 

3.1.2 Supersonic Formulation 

To develop a supersonic quadrupole formula- 
tion, equation (15) is transformed from a collapsing 
sphere formulation to an emission-surface formula- 
tion. The emission-surface formulation is appropri- 
ate for supersonic-source motion because it does not 
have a Doppler singularity |1 — M r | in the denomina- 
tor of the integrand. Farassat and Brentner 14 took 
the time derivatives inside the integrals by recogniz- 
ing that the entire rotor plane can be considered as 
the emission-surface, hence the limits of integration 
are not a problem. Next they changed coordinates 
from a frame fixed to the undisturbed medium to 
a frame always aligned with the rotor blade, (i.e. 
(x,£) — y (r) ,t)). The effect of this operation on the 
tensor Q with components Qij is 


dm 

dt 


ret 


= [^rQ«j] ret 


ret 


(17 


where 77 is the position vector in the rotating frame 
and r is the source time. Here V = dr] /dr is the 
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velocity of the point 77 specified in the frame fixed 
to the undisturbed medium. We note that V has 
no component normal to the rotor plane. It is im- 
portant to recognize that when we refer to Qij \ v we 
really mean that the components of the tensor 
represented in coordinates that are instantaneously 
aligned with the rotating frame. Thus equation (17) 
provides the time derivative of Qij\ x in the station- 
ary frame in terms of Qij\ n which is specified in the 
coordinates of the moving frame. Using the operator 
notation L T , the final emission-surface formulation 
may be written 

4 7T P ' Q (x,t) = ^ J \hfi[ L lQij}ret d S 

+ >/ ^2 ffiiTjLrQij ~ L T Qu] ret dS 

+ J ^3 [SQrr ~ Qii] ret • (18) 

Notice that the operator L T only operates on 
because r* and fj do not depend upon t or r. We 
refer to equation (18) as formulation Q2. Examining 
the definition (17) and equation (18), we realize that 
formulation Q2 is a singularity free expression for su- 
personic quadrupole noise. Note that equation (18) 
has second space and time derivatives of Qij as well 
as first space derivatives in the rotor plane. These 
quantities are available in the CFD postprocessor 
that is used to compute Qij for acoustic calculations. 
As it stands, formulation Q2 is valid for subsonic 
and supersonic quadrupole noise prediction for heli- 
copter rotors in hover or forward flight. (Full details 
are given in reference 14.) 

Formulation Q2 has been implemented in a 
demonstration code known as WOPWOP2+ . WOP- 
WOP2+ differs significantly from WOPWOP+ 9-11 
in that it uses an emission-surface formulation to 
compute thickness and loading noise, as well as the 
quadrupole noise. The construction of the emission 
surface and subsequent integration over the emission 
surface is performed using the method of march- 
ing cubes integration developed by Brentner. 15 In 
reference 14, it is shown that the subsonic (WOP- 
WOP+) and the supersonic (WOPWOP2+ ) formu- 
lations give identical results when the quadrupole 
grid extent is the same. Ianniello 16 has also devel- 
oped quadrupole noise prediction codes which uti- 
lize the far-field approach of Brentner 9, 10 and he has 
developed a sophisticated emission-surface construc- 
tion and integration scheme. 17 

We now present HSI noise calculations for a two- 
blade model-scale UH-1H rotor tested in hover with 
tip Mach numbers 0.88, 0.9, 0.925, and 0.95 (See 
reference 18 for test information). The quadrupole 


grid in these computations extends 1.867? beyond 
the blade tip for all the WOPWOP2+ calcula- 
tions shown in figure 3. For comparison, we have 
also shown the signature predicted by WOP WOP + 
which includes quadrupole sources only up to the 
sonic circle. The agreement of the WOPWOP2-h 
signature with the measured data is excellent and 
better than that of WOPWOP+ for each case. For 
the more intense cases ( Mh > 0.90), the agreement 
of the WOPWOP2+ prediction with the measured 
acoustic pressure signature is not fully satisfactory 
because the WOPWOP2+ prediction overpredicts 
the negative peak pressure. This is apparent in fig- 
ure 3 for the times between the WOP WOP + and the 
WOPWOP2+ shock locations. Farassat and Brent- 
ner 14 identify this over predict ion as numerical oscil- 
lation related to the bifurcation of the quadrupole 
source region at a critical source time. 

3.2 FW-H equation revisited 

One example of developing an alternative but 
equivalent source description has been the utilization 
of the FW-H equation on a permeable surface. In 
rotor noise prediction, the surface / = 0 has usually 
been assumed to be coincident with the rotor blade 
surface and impenetrable ( u n = u n ). A relaxation of 
that assumption is useful because it allow consider- 
ation of either conveniently placed fictitious surfaces 
or physical surfaces which permit flow through them. 

For a permeable surface (fictitious or physical), 
equation (6) is the appropriate expression of the 
FW-H equation. It is clear from both the origi- 
nal Ffowcs Williams and Hawkings paper 1 and later 
work [19, Chapter 11, Sec 10] that Ffowcs Williams 
understood the value of utilizing equation (6) on a 
permeable surface. Recently di Francescantonio im- 
plemented this form of the FW-H equation for rotor 
noise prediction. 20 At nearly the same time, Brent- 
ner and Farassat demonstrated the relationship be- 
tween the FW-H equation and the Kirchhoff equa- 
tion for moving surfaces. 21 Many others have fol- 
lowed suit by quickly adopting the utilization of the 
FW-H equation on a permeable surface. 22-24 

One question to be considered is “What is the 
advantage of a fictitious surface that is not coinci- 
dent with the physical body?” Although it may not 
be immediately obvious from equations (6) or (7), 
any physical acoustic sources enclosed by the / = 0 
surface only contribute through the surface-source 
terms. Any physical sources of sound or propaga- 
tion effects outside / = 0 only contribute through 
the volume source term. Hence, if we can enclose all 
physical sources inside / = 0 we have no contribu- 
tion from the volume source — the quadrupole can be 
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t, msec t, msec 

Figure 3. Comparison of WOPWOP+ ( ) and WOPWOP2+ ( ) predicted acoustic pressure with ex- 

perimental data 18 (□) for hovering model UH-1H rotor. Quadrupole grid in WOPWOP+ prediction extended 
almost to sonic circle and in WOPWOP2-f predictions extended 1.86i? beyond the rotor tip. 


legitimately neglected. Without the quadrupole, a 
significant computational savings is realized because 
the volume integration is no longer required. Fur- 
thermore, the amount of flow-field data required for 
a surface integration is much less than for an inte- 
gration of the the volume surrounded by the surface. 

The enabling key to utilization of the FW-H equa- 
tion on a permeable surface is the availability of an 
accurate flow-field description from CFD. Only re- 
cently has CFD matured point where it can provide 
sufficiently accurate, unsteady flow-field data on the 
integration surface. Even so, CFD computations for 
a rotor in forward flight are extremely demanding. 
Thus the coupling of the FW-H equation and CFD 
provides a mutually beneficial approach to comput- 
ing the noise. The CFD computation only is needed 
in the acoustic source region — not all the way to the 
observer — and the FW-H equation provides an effi- 
cient method of predicting the sound field away from 
the source region. Through the utilization of the 


permeable surface formulation of the FW-H equa- 
tion, the acoustic calculations can be made compu- 
tationally efficient even for complicated, nonlinear 
acoustic sources. 

To demonstrate the robustness and accuracy of 
the permeable surface application of the FW-H 
equation, a comparison between the predicted acous- 
tic pressure and experiment for a hovering rotor is 
shown in figure 4. In the figure, the result for a 
WOPWOP+ prediction is also shown. The same 
CFD calculation was utilized for the input data all 
the predictions. Two permeable surface FW-H com- 
putations are show in figure 4: first, an integration 
surface coincident with the rotor blade to predict 
thickness and loading noise; and second, an integra- 
tion surface located approximately 1.5 chordlengths 
away from the blade to predict the total noise. Note 
that the thickness noise predictions from WOP- 
WOP+ and FW-H are identical and there is only 
a small difference in the predicted loading noise. 
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time, msec 


Figure 4. Comparison of noise components predicted 
by the FW-H/RKIR and WOPWOP+ codes for a 
hover UH-1H model rotor (Mh = 0.88, inplane ob- 
server 3.4i? from rotor hub). 


The total noise, which includes the effect of the 
quadrupole, is also in very close agreement even 
though the volume used in WOP WOP + is not iden- 
tical to the region enclosed in the FW-H permeable- 
surface integration. 


3.3 Problems with the Klrchhoff formulation 

The Kirchhoff formulation for moving surfaces is 
an alternative formulation that has been used for ro- 
tor noise prediction over the past decade. The Kirch- 
hoff formulation gained rapid acceptance shortly af- 
ter its publication by Farassat and Myers in 1988 25 
precisely because it appeared to offer the same ben- 
efits just presented for the permeable-surface appli- 
cation of the FW-H equation. Nevertheless, recent 
work by Brentner and Farassat 21 and Singer et al. 26 
have shown that the Kirchhoff formulation is unre- 
liable for aeroacoustic problems in practice. This 
observation is sufficiently important that some ex- 
planation is desirable. 

An embedding procedure similar to that used to 
derive the FW-H equation above was applied to 
the wave equation by Farassat and Myers to derive 
the Kirchhoff formulation for moving surfaces. The 
generalized wave equation — which is the governing 


equation for the Kirchhoff formulation — becomes 


DV(M) . 


= Qki 


(19) 


where M n = v n /c. In this equation p ! must be com- 
patible with the wave equation, hence, equation (19) 
is valid only in the region of the fluid in which the 
wave equation is the appropriate governing equation. 

Through the utilization of the continuity and mo- 
mentum equations we can rewrite the permeable sur- 
face form of the FW-H equation (6), as 

□ V(x,f) = Qkir + fc-dl 

+ - C 2 p']^-S(f) + | [V - c*p')^6(f)\ 

- 7^7 [pUiUj]niS{f) - -j^[puiU n 6(f)\ . (20) 

This form of the FW-H equation highlights the 
differences between the Kirchhoff formulation and 
the FW-H equation. Notice that all the additional 
source terms in equation (20) are second-order terms 
in perturbation quantities outside of the source re- 
gion (i.e. p l — c 2 p f -C 1 and pUiUj 1). For lin- 
ear wave propagation, each of these terms would be 
identically zero and the Kirchhoff formulation and 
the FW-H equation would be in complete agree- 
ment. (A more detailed comparison and discussion 
are found in reference 21.) 

To illustrate how the Kirchhoff formulation can 
give misleading results, consider a series of cylin- 
drical integration surfaces which enclose a hovering 
rotor blade, as shown in figure 5. In figure 6, the 
acoustic pressure for an in-plane observer has been 
computed using both the FW-H equation and the 
Kirchhoff formulation for each of the integration sur- 
faces. 

Both methods agree reasonably well with the data 
when the integration surfaces are more than about 
0.7 chords from the blade surface. The predicted 
acoustic pressure from the Kirchhoff computations 
for integration surfaces that are closer to the blade 
are unrealistic. The acoustic pressure predicted by 
the FW-H equation, however, is well behaved and 
is modified only by the fact that not all of the 
“quadrupole” noise source is included when the inte- 
gration surface is too close to the blade and the vol- 
ume source term in equation (6) has been neglected. 

A second example illustrates another problem that 
can occur with the Kirchhoff formulation. In this 
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Figure 5. Concentric cylindrical integration surfaces 
used for noise computation for a hovering rotor. 



k=2 

k=7 

k=12 

k=18 

k=21 

data 


(a) 



Figure 6. Predicted acoustic pressure for various 
integration surface locations for an observer located 
SAR from a UH-1H model rotor hovering at Mh = 
0.88. The experimental data (□) is from reference 
18. (a) FW-H prediction; (b) Kirchhoff prediction 



Figure 7. Vorticity field computed from CFD. 
FW-H integration surfaces are atr = 0.5D,r = 
1.5D, r = 2.5Z>, and r = 5 . ID 

Kirchhoff FW-H 



t, sec t, sec 


Figure 8. Acoustic signals computed for various in- 
tegration surfaces that correspond to those indicated 
in Figure 7. 


case, consider the viscous flow past a circular cylin- 
der shown in figure 7. In this situation it is ex- 
pected that the vortices shed by cylinder have a very 
small contribution to the sound produced, hence the 
acoustic signal should be relatively unaffected by the 
placement of the integration surface. Figure 8 shows 
that this is indeed the case for the FW-H computa- 
tion, but the Kirchhoff computation does not even 
converge to a value but is entirely erroneous. This 
is significant for rotor noise prediction because rotor 
wakes inevitably must pass through the integration 
surface. 

These two examples demonstrate numerically that 
the Kirchhoff formulation is not reliable for rotor 
noise prediction. Fortunately, from a computational 
point of view, there is very little difference between 
the methods. All of the computational advantages 
originally sought from the Kirchhoff method are 
available using the FW-H equation. Furthermore, 
all of the physical insight that the FW-H equation 
traditionally has provided is also still available. 
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4. Summary 

In the last several years significant progress has 
been made in developing efficient quadrupole noise 
prediction formulations. Much of the work was 
motivated by the need for an efficient high-speed 
impulsive-noise prediction tool. The FW-H equa- 
tion has proven to be a fruitful source of integral 
formulations and the desire to avoid the compu- 
tation of the quadrupole noise altogether has lead 
to the utilization of the FW-H equation on perme- 
able surfaces. The permeable-surface FW-H equa- 
tion actually embodies all of HSI prediction meth- 
ods because the full formulation still includes the 
quadrupole outside of the integration surface. It has 
also been demonstrated that the Kirchhoff formu- 
lation, while fine for acoustic problems, can be un- 
reliable for aeroacoustics. Fortunately the FW-H 
equation is equally efficient and much more robust. 

Although this paper has not dealt directly with 
rotorcraft aerodynamic computations, it should be 
pointed out that the acoustic formulations dis- 
cussed require highly accurate solutions — both spa- 
tially and temporally — as input data. The aerody- 
namic calculations required for HSI noise in forward 
flight are extremely challenging. An accurate, first- 
principles, blade- vortex interaction computation is 
still beyond the ability of present day computers and 
CFD. Nevertheless, progress is being made and the 
acoustic propagation theory and codes will be up to 
the task. 
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